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Abstract 

We present a preliminary determination of the inclusive branching fraction for the rare radiative 
penguin transition B — > X s ^. The measurement is based on a data sample of 60 million BB 
pairs collected between 1999 and 2001 with the BABAR detector at the PEP-II asymmetric-energy 
e + e~ B Factory at SLAC. We study events containing a high-energy photon from one B (or 
B) decay and a tagging primary lepton from the decay of the other B meson. By this means, 
we are able to reduce a significant component of the background without introduction of model 
dependent uncertainties in the event selection efficiency. We determine the branching fraction 
BiB -> X s7 ) = [3.88 ± 0.m(stat.) ± 0.37(syst.) ±gj| (model)} x 10" 4 , which is consistent with 
Standard Model predictions and provides a constraint on possible new physics contributions to the 
electromagnetic penguin amplitude in B decays. 
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1 Introduction 



The parton-level b — > S7 "radiative penguin" transition rate can be calculated at next-to-leading 
order to a precision of 10% in the Standard Model (SM) . The presence of non-SM particles in the 
virtual loop mediating this transition may lead to substantial deviations from this predicted rate 
and has been the subject of numerous theoretical investigations Q. The inclusive B — > X s *y rate is 
equal to the calculated parton-level b —>■ rate according to quark-hadron duality [|J . There have 
been several measurements of the branching fraction for B — ► X s ^/ to date PJ. The confinement of 
the b quark in the B meson results in model dependent assumptions for the E~ spectrum, while the 
X s system fragments non-perturbatively into the detectable final state particles. Thus, in order 
to access the parton model process, it is desirable to impose as few requirements as possible on 
the Xg'j system, motivating a fully-inclusive approach. This, however, makes the measurement 
susceptible to large backgrounds, which must be suppressed without making such requirements. 

In this paper we present a measurement of the branching fraction B(B — » X s ^) in which a 
significant component of the background is removed by making requirements on the other B meson 
in the event, rather than on the signal B — * X s -f decay. The remaining background is constrained 
by independent control samples. Although a requirement on is still necessary, the accessible 
energy range is limited by the statistics of the control samples, which scale directly with the signal 
size. Hence this technique scales well with the increasing data sample anticipated at BABAR. This 
fully-inclusive technique does not distinguish between decays of charged and neutral B mesons, so 
both are used. For simplicity we refer to both B and B as "B mesons". 

2 The BABAR detector and dataset 

The data were collected with the BABAR detector at the PEP-II asymmetric-energy e + (3.1 GeV) - 
e~(9GeV) storage ring. A description of the BABAR detector can be found in Ref. Q. Charged par- 
ticles are detected and their momenta measured by a combination of a silicon vertex tracker (SVT), 
consisting of five double-sided layers, and a central drift chamber (DCH), in a 1.5-T solenoidal field. 
We identify leptons and hadrons with measurements from all detector systems, including the energy 
loss (dE/dx) in the DCH and SVT. Electrons and photons are identified by a Csl electromagnetic 
calorimeter (EMC). Muons are identified in the instrumented flux return (IFR). A Cherenkov ring 
imaging detector (DIRC) covering the central region provides particle identification. 

We use Monte Carlo simulations of the BABAR detector based on GEANT 4.0 || to optimize our 
selection criteria, to determine signal efficiencies and to estimate part of the background component. 
Events taken from random triggers are mixed with simulated events to take beam backgrounds and 
some varying detector conditions into account. 

The results in this paper are based upon an integrated luminosity of 54.6 ± 0.8 fb -1 of data, 
corresponding to (59.6±0.7) x 10 6 BB meson pairs recorded at the T(AS) resonance ( "on-resonance" ) 
and 6.40 ± 0.10 fb -1 at 40MeV below this energy ("off-resonance"). The off-resonance data, which 
constitute a fraction f g = 10.5% of the total data set, are used for background estimation. The 
number of BB meson pairs is determined from the excess hadronic events relative to muon pairs 
in on-resonance data compared to off-resonance data [||]. 
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3 Analysis method 



The B — > X s j decay signature is characterized by the high energy photon. The background to 
the B — > X s 7 signal consists of two components. First, it can arise from other B meson decays, 
in which the photon candidates are predominantly from 7r° or i] (or rarely lo) decays, or - in a 
small fraction of cases - from a hadronic interaction of a neutron or K® in the electromagnetic 
calorimeter. Second, it can arise from continuum qq or r + r~ production, where q can be a u,d,s or 
c quark, with the high-energy photon originating either from initial-state radiation (ISR) or from 
similar sources to those in B decays. Figure [l] demonstrates the scale of the background rejection 
problem (based on Monte Carlo simulation). 




1.8 2 2.2 2.4 2.6 2.8 3 3.2 3.4 

E Y (GeV) 

Figure 1: The energy distribution, in the T(4S) center of mass, of simulated photon candidates, 
after "photon quality" cuts designed to reduce backgrounds from 7r°s, 77s and hadrons, but before 
other selection requirements. Shown are B — > X s ^f signal (dark shading), BB background (grey 
shading) and continuum background (unshaded), all normalized to 54.6 fb -1 . The high-end tail for 
BB is mainly due to residual hadrons. 

The method for extracting the signal from data is to subtract the continuum background based 
on off-resonance data, and the BB contribution based on Monte Carlo predictions, where the 
latter are validated (and if necessary corrected) using our control samples. A series of selection 
requirements is made to greatly suppress backgrounds (primarily continuum) before the subtraction. 
The selection criteria for this analysis are optimized for statistical precision by maximizing ns 2 /(ns+ 
ne + nc//off). The number of signal candidates ng is estimated using a Monte Carlo simulation 
that incorporates the model of Kagan and Neubert [ffl] for the spectrum and assumes a central 
theoretical value of B{B — > X s ~/) = 3.45 x J|]. The BB and continuum backgrounds, ne 

and nc, are estimated with Monte Carlo simulation. The signal and background expectations are 
normalized to 54.6 fb _1 . The continuum is additionally weighted by l// ff in the optimization to 
allow for the smaller size of the off-resonance data set. 

Initially we require a high-energy photon candidate with 1.5 < E* < 3.5 GeV in the e + e~ 
center-of-mass frame. (An asterisk denotes quantities computed in the e + e~ center-of-mass frame, 
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as opposed to the laboratory frame.) This is a loose selection that is tightened later to define the 
signal region and also control "sideband" regions. A photon candidate is defined as a localized 
energy maximum Q in the calorimeter acceptance —0.74 < cos# < 0.93, where 9 is the lab-frame 
polar angle relative to the e~ beam direction. It must be isolated by 40 cm at the inner surface 
of the EMC from any other photon candidate or track. We veto candidates whose lateral energy 
profile is inconsistent with a single photon shower; these can arise from hadronic interactions in the 
calorimeter, or from 7r°s in which the decay photons have merged into one shower ("merged 7r°s"). 
In addition we veto photons from a vr°(r/) when the invariant mass of the combination with any 
other photon of energy greater than 50(250) MeV is within ~ 2.7(2.2) sigma of the nominal vr°(?7) 
mass, 115(508) < M 77 < 155(588) MeV/ c 2 . However, we retain such vetoed photons as a control 
sample to validate the Monte Carlo modeling of the BB background. The photon candidates that 
fail the lateral profile cut form the "hadron anti-veto" sample (it also includes a small component 
of merged vr°s); while the photons that pass the lateral profile cut, but are rejected by the 7r° or 
rj veto, form the "7r°(ry) anti-veto sample". All three independent samples, the signal and the two 
anti-veto samples, are required to pass all the subsequent selection criteria. 

To suppress continuum backgrounds we require a "lepton tag". This is a high-momentum 
electron (p* > 1.3GeV/c) or muon (p* > 1.55GeV/c). For B — > X s ^ signal events the lepton arises 
from the semileptonic decay of the other B meson. Leptons also occur in the continuum background, 
most notably from the semi-leptonic decays of charm hadrons, but are produced significantly less 
frequently and with lower momentum than from a B decay. An electron candidate must have 
a ratio of calorimeter energy to track momentum, an EMC cluster shape, a DCH dB/dx and a 
DIRC Cherenkov angle (if available) consistent with an electron. A muon candidate must satisfy 
requirements on the measured and expected number of interaction lengths penetrated, the position 
match between the extrapolated DCH track and IFR hits, and the average and spread of the number 
of IFR hits per layer. The lepton candidates that do arise from continuum events are often pions 
misidentified as muons or, less frequently, electrons. Fake leptons are usually not associated with 
an undetected neutrino, unlike real leptons from a semi-leptonic B decay. We therefore require 
that the missing energy in the event be greater than 1.2 GeV. Further discrimination is obtained 
by considering the angle between the lepton and the high-energy photon. For signal events, the 
cosine of this angle, cos#* e (cos# 7At ), has an approximately flat distribution, except for a small 
peak at cos 9*^ = — 1 arising from cases where the muon candidate and photon both originate 
from the signal B decay. (This can occur when a pion from the X s decay fakes a muon.) In 
contrast, the continuum background is strongly peaked at cos #* e (cos 6* ) = — 1 as a consequence 
of the jet-like topology, with the peak at cos #* e (cos 9*) = +1 suppressed by the photon isolation 
requirement. We require cos 9* e (cos 9*) > —0.75 (—0.70). This also gives a small suppression 
of the BB background for cases where the high energy photon and lepton arise from the same 
B decay. Note that, since the tagging requirement is imposed on the "other" B rather than on 
the signal B, we are able to reject continuum background without imposing requirements on the 
signal decay. The signal efficiency of the combined lepton tag, missing energy and cos #* e (cos 9*^) 
requirements is 5%, while the continuum background is suppressed by a factor of approximately 
1200. The momentum-dependent electron and muon efficiencies used for these and other Monte 
Carlo results have been measured in data with a sample of radiative Bhabha events and samples 
of e + e~ — > fi + fi~j and e + e~ — » e + e~^ + /i~ events, respectively. 

Further continuum rejection is gained by imposing selection criteria on the event topology. In 
the center-of-mass system BB decays are isotropic, while the inclusive ir (rj) and ISR components 
of the continuum background have a two and three-jet topology, respectively. We compute the 
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ratio of the second to the zeroth order Fox- Wolfram moment || in both the center-of-mass frame, 
i?2, and the rest frame of the system recoiling against the high energy photon, R 2 . The recoil frame 
recovers the two-jet topology for ISR events. We require R 2 < 0.45 and R 2 /R 2 < 1- We compute 
the energy sums of all particles (except the photon candidate) whose momentum vectors lie within 
0° — 30° and 140° — 180° cones about the photon direction, referred to as E% and E£, respectively. 
We require E) < 1.1 GeV and 1.6 < E* h < 3.6 GeV. 

Figure || shows the Monte Carlo prediction for the E* distribution for signal, BB background 
and continuum background after the preceding selection criteria, and the expected E* signal dis- 
tribution for different model parameters. The signal region is defined to be 2.1 < E* < 2.7 GeV, 
while the regions 1.7 < E* < 1.9 and 2.9 < E* < 3.5 GeV, which are dominated by BB and 
continuum background, respectively, serve as "control" regions for the estimation of these back- 
ground components. The choice of the final E* signal region, in particular the lower bound on 
E* is a balance between the systematic uncertainty of the estimated BB component, which is the 
dominant systematic uncertainty of the measurement, and the model dependence. As the lower 
E* bound is reduced, the model dependence decreases but the BB background, which is estimated 
with a Monte Carlo simulation, rises sharply. 




E y (GeV) E* Y (GeV) 

Figure 2: Left: The reconstructed E* distribution expected from Monte Carlo simulation after 
the selection criteria. The B — > X s j signal assuming B(B -> X s -y) = 3.45 x 1(T 4 (dark shading), 
BB background (grey shading) and continuum background (unshaded) are normalized to 54.6 fb _1 . 
Right: The generated E* spectrum before cuts (arbitrary normalization) for different values of the 
b quark mass mi,, using the model of Kagan and Neubert |ffl|. Our signal region is defined for the 
corresponding reconstructed quantity as 2.1 < E* < 2.7 GeV. 



4 Systematic studies 

In Table |l] we list the fractional systematic uncertainties on the extracted branching fraction B(B —> 
X s j). The dominant uncertainty is the modeling of the BB background. To test the modeling of 
the E* spectrum of the BB background component, we compare the two anti-veto samples, after 
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Table 1: The systematic error expressed as a percentage of B(B — > X s j). The total is the quadratic 
sum of all contributions, 



Systematic Uncertainty 


AB(B -> X 8 'YVBf%) 


BB background estimation 


1 ^7 Q 


Signal Monte Carlo statistics 


±Z.D 


Photon efficiency 


4-9 ^ 
±Z.O 


Lepton tag efficiency 


4-9 n 




4-1 n 




+1 n 

_1_ ± ,\J 


Photon distance requirement 


±2.0 


7r°/?7 veto 


±1.0 


BB count 


±1.1 


Total 


±9.3 


Signal model dependence 


+11.2, -5.9 



the continuum contribution has been subtracted using the off-resonance data, with generic BB 
Monte Carlo event samples that have passed the same selection criteria described above. Figure |3| 
shows the comparison between Monte Carlo simulated events and data for the 'V and r\ anti-veto 
sample" and the "hadron anti-veto" sample. To take any differences into account, we separately fit 
the ratio of the measured data to Monte Carlo prediction for each control sample with a first-order 
polynomial, and compute the average of the fit over the range 2.1 < E* < 2.7 GeV, weighted by the 
E* spectrum expected for the corresponding background to our signal sample. We then correct the 
E* spectrum of the BB Monte Carlo simulation used to estimate the background in the signal region 
according to these fit results, further weighted by the relative contributions of the 7r° and rj (91%) 
and hadronic and merged ir° components (9%) to the background. To estimate the uncertainty in 
this correction we find the variation incurred when using a zero-order polynomial fit rather than 
first-order, or moving the lower bound of the fit region from 1.7 to 1.8 GeV, or simply considering 
the ratio of data to Monte Carlo events in the region 2.1 < E* < 2.7 GeV. This variation is summed 
in quadrature to the statistical error of the first-order polynomial fit. The correction factor for the 
total BB background in the region 2.1 < E* < 2.7 GeV is 0.89 + 0.16. The final uncertainty on the 
BB background estimate is computed by counting the expected number of Monte Carlo events and 
then correcting as described, with the uncertainty being the sum in quadrature of the statistical 
error of the Monte Carlo expectation and the uncertainty in the correction factor. This results in 
a corrected expectation for BB events in 2.1 < E* < 2.7 GeV of 222 ± 43. (This is translated to 
the fractional entry in Table |l| by dividing the uncertainty by the extracted signal from Section 0.) 

The remaining systematic uncertainties are from the uncertainty in the signal efficiency due 
to the photon selection and the lepton tag requirements, plus a small contribution from BB 
event counting. The photon efficiency is measured by comparing the ratio of events Nfr^ — > 
h^ z TT°)/N(T :iz — ► h^ir !? ) to the previously measured branching fractions Q. The photon isolation 
and 7r°/?7 veto efficiency are dependent on the event multiplicity. Their simulation is tested by 
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Figure 3: Left: The E* distribution for the 7r° and 77 anti-veto control sample for data (points) 
compared to the Monte Carlo prediction (solid line). Right: Corresponding distribution for the 
hadron anti-veto control sample. 



"embedding" Monte Carlo-generated photons into both an exclusively reconstructed B meson data 
sample and a generic B meson Monte Carlo sample. The photon-energy resolution is measured in 
data using ir° and rj meson decays and a sample of virtual Compton scattering events. The energy 
scale uncertainty is estimated with a sample of 77 meson decays with approximately equal-energy 
photons; the deviation in the reconstructed 77 mass from the nominal 77 mass provides an estimate of 
the uncertainty in the measured single photon energy. Lepton tag systematics have been estimated 
by coherently varying the measured efficiencies by ±lcr. The total systematic uncertainty of 9.3% 
is the quadratic sum of all the contributions itemized in Table |]. 

We use an E* spectrum based on the model of Kagan and Neubert Q to determine the signal 
efficiency. In this model, the E^ spectrum, which is dual to the mass spectrum of the X s system 

9 9 
fry, — 771 

through the relation Ej = B 2 mB Xs , where E^ is the energy of the photon in the B meson rest frame, 
has two components. The region mj s < m cu t ff is described by a relativistic Breit-Wigner for the 
B — > K* (892)7 decay. The region mx a > Tn cuto g is described by a spectrum parameterized in terms 
of the 6-quark mass mb- The efficiency is sensitive to the value of m& because that determines how 
much of the spectrum falls outside 2.1 < E* < 2.7 GeV (right plot in Figure |||), but insensitive to 
the fraction of if* (892) or to the m cu t fr value. Kagan and Neubert pj recommend varying m& from 
4.65 to 4.95 GeV/ c 2 to estimate model dependent variations. We therefore compute the efficiency 
of the signal using mj, = 4.80 ± 0.15 GeV/c 2 . We fix m cuto g- to 1.1 GeV/c 2 , the value measured in a 
BABAR semi-inclusive analysis of B — ► X s ^f ||, and set the if* (892) fraction equal to the integral 
of the discarded continuum spectrum below m cu t ff for each m&, a procedure suggested by Kagan 
and Neubert, to find resulting shifts of +11.2% and -5.9% in B{B — > Xg'y). To test the insensitivity 
of the result to the details of the relative contributions of B — > K* (892)7 and X s assumed in the 
model, we explicitly vary m cu t ff by ±100MeV/c 2 , and independently the fraction of B — ► if* (892)7 
by a factor of two. We find negligible additional uncertainty on B(B — > X s j). 
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5 Results 



In order to reduce experimenter bias, the region 1.9 < E* < 2.9 GeV was kept hidden until all 
selection criteria had been finalized and backgrounds estimated. Figure || shows the E* spectrum 
for data compared to the total predicted background. The control regions 1.7 < E* < 1.9 GeV and 
2.9 < E* < 3.5 GeV show no statistically significant excess of data, as expected, but a signal is 
plainly visible in between. In the signal region, 2.1 < E* < 2.7 GeV, we find 543 net signal events. 
Note that the selection procedure does not discriminate against the expected small component of 
b — > dj, so these events are included in our total. Combined with a selection efficiency within 
the E* range of 1.28 ± 0.06% and the dataset size of (59.6 ± 0.7) x 10 6 BB meson pairs, this 
yield corresponds to a measurement of the partial branching fraction for 2.1 < E* < 2.7 GeV of 
B{B -> Xi) = [3.55 ± 0.32(atut.) ± 0.32(syst.)] x 10" 4 . For a Kagan and Neubert-based model 
with m& = 4.80 GeV/ c 2 (see Fig. ||), the overall efficiency (including the effect of the E* cuts) is 
1.127 ± 0.055%, resulting in B(B -» Xi) = [4.05 ± 0.37(stat.) ± 0.38(syst.) ±%f 4 {model)] x 10" 4 . 
Finally, we subtract from this a contribution from b — > dj. In the Standard Model the theoretical 
expectation is B{B — > X c ij)/B(B — ► Xgj) = \V t d/V ts \ 2 . Assuming the signal efficiency for B — > X^ 
to be equal to that for B — > X s j and |14rf/14 s | = 0.20±0.04 we then scale the measured branching 
fraction by 0.96 ± 0.02 to give the preliminary result: 

B(B -» X sl ) = [3.88 ± 0.36(stat.) ± 0.37{syst.) ±g;|| (model)] x 10~ 4 . 



900 




E y (GeV) 

Figure 4: The E* distribution of on-resonance data (solid points) compared to background expec- 
tation. All errors are statistical only (including, just for this figure, BB background statistics). 
Bin-by-bin systematic uncertainties and correlations have not yet been studied; the systematics 
quoted in the text apply only to an integral measurement from 2.1 to 2.7 GeV. 
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6 Summary 



We find a preliminary branching fraction for B — > X s ^ of [3.88 ± 0.36(stat.) ± 0.37(syst.) ±023 
(mode/)] x 10 -4 using an inclusive technique. Figure |5| compares this result with theoretical pre- 
dictions and with previous measurements. We do not see any evidence for a departure from the 
next-to-leading order Standard Model predictions. The systematic precision is limited by the size 
of the BB background control samples, which scale in proportion to the signal sample. The 
systematic precision in turn limits the lower bound, E* > 2.1 GeV, of the signal region. As the 
larger datasets anticipated at the B-factory become available, the systematic uncertainty in the 
BB background will be reduced along with statistical uncertainties. Furthermore, this reduction 
in the systematic uncertainy may allow using a lower minimum-B* bound for the signal region, 
which will in turn lead to smaller model dependence. 
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Figure 5: The BABAR measurement compared to previous experiments || and to theoretical pre- 
dictions Iffll. 
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